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Abstract. We study the dependence of the cooling of isolated neutron stars on the magnitude of the 3P2 neutron 
gap. It is demonstrated that our "nuclear medium cooling" scenario is in favor of a suppressed value of the 3P2 
neutron gap. 
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1. Introduction 

As the result of many works the so called "standard" 
scenario of neutron star (NS) cooling emerged, where 
the main process responsible for the cooling is the mod- 
ified Urea process (MU) nn — > npev calculated us- 
ing the free one pion exchange between nucleons, see 
|Friman fc Maxwell (1979) This scenario explains only the 
group of slow cooling data. To explain the group of 
rapid cooling data the "standard" scenario was supple- 
mented by one of the so called "exotic" processes either 
with pion condensate, or with kaon condensate, or with 
hyperons, or involving the direct Urea (DU) reactions, 
see |Tsuruta (1979)] |Shapiro fc Teukolsky (1983)| and ref- 
erences therein. All these processes may occur only for 
densities higher than a critical density, (2 -j- 6) n , de- 
pending on the model, where no is the nuclear saturation 
density. 

Then the pair breaking and formation (PBF) processes 
permitted in nucleon superfluids have been suggested. 



and 3P2 superfluids. Schaab et al. (1997) have shown that 
the inclusion of the PBF processes into the cooling code 
may allow to describe the " intermediate cooling" group of 
data (even if one artificially suppressed medium effects). 
Thus the "intermediate cooling" scenario arose. Then the 
PBF processes were incorporated in the cooling codes of 
other groups, elaborating the "standard plus exotics" sce- 
nario, see |Tsuruta et al. (2002)[ [Yakovlev et al. (2004)[ 
Page et al. (2004)| Some papers included the possibility 
of internal heating that results in a slowing down of 
the cooling of old pulsars, see Tsuruta (2004) and Refs. 
therein. However, paying the price for a simplification of 
the consideration, calculations being performed within the 
"standard plus exotics" scenario, did not incorporate in- 



Flowers et al. (1976) calculated the emissivity of the ISq 



neutron pair breaking and formation (nPBF) process 
and Voskresensky & Senatorov (1987) considered a gen- 
eral case. Neutron and proton (pPBF) pair breaking and 
formation processes were incorporated within a closed dia- 
gram technique including correlation effects. Up to factors 
of the order of 1 numerical estimates are valid both for ISq 



medium effects. Recently Page et al. (2004) called this ap- 
proach the "minimal cooling" paradigm. 

The necessity to include in-medium effects into the NS 
cooling problem is a rather obvious issue. It is based on 
the whole experience of condensed matter physics, of the 
physics of the atomic nucleus and it is called for by the 
heavy ion collision experiments, see |Migdal et al. (1990)1 
Rapp fc Wambach (1994)[|Ivanov et al. (2001)| The rele- 
vance of in-medium effects for the NS cooling problem has 



been shown by Voskresensky & Senatorov (1984), (1986), 
(1987), |Migdal et al. (1990)||Voskresensky (2001)| and the 
efficiency of the developed "nuclear medium cooling" sce- 
nario for the description of the NS cooling was demon- 
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and then by Blaschke et al. (2004) In that paper, it was 



shown that it is possible to fit the whole set of cooling 
data available by today. Besides the incorporation of in- 
medium effects into the pion propagator and the vertices, 
it was also exploited that the 3P 2 neutron gaps are dra- 
matically suppressed. The latter assumption was moti- 
vated by the analysis of the data (see Figs 12, 15, 20 - 
23 of Blaschke et al. (2004) I and by recent calculations of 



the 3P2 neutron gaps by Schwenk & Friman (2004) 




n / n. 



Fig. 1. Neutron and proton pairing gaps according to 
model I (thick solid, dashed and dotted lines) and ac- 
cording to model II (thin lines), see text. The ISq 
neutron gap is the same in both models, taken from 
Ainsworth et al. (1989) 



I and II within the "nuclear medium cooling" scenario. 
As it was checked there, since the ISo neutron pairing 
gap exists only within the crust, dying for baryon den- 
sities n > 0.6 no, its effect on the cooling is rather mi- 
nor. Opposite, the effect on the cooling arising from the 
proton ISo pairing and from the neutron 3P 2 pairing, 
with gaps reaching up to rather high densities, is pro- 
nounced. The NS cooling essentially depends on the val- 
ues of the gaps and on their density dependence. Findings 
of |Schulze et al. (1996)||Lombardo fc Schulze (2000)[ who 



incorporated in-medium effects, motivated us to check the 
possibility of rather suppressed 15*o neutron and proton 
gaps. For that aim we introduced pre-factors for 15*0 neu- 
tron and proton gaps which we varied in the range 0.2 1, 



see Figs. 18 and 19 of Blaschke et al. (2004) 




6.4 6.6 
log 10 (T[K]) 



In spite of many calculations which have been per- 
formed, the values of nucleon gaps in dense NS matter 
are poorly known. This is the consequence of the expo- 
nential dependence of the gaps on the density depen- 
dent potential of the in-medium NN interaction. This 
potential is not sufficiently well known. Gaps that we 
have adopted in the framework of the "nuclear medium 
cooling" scenario, see Blaschke et al. (2004) are presented 



in Fig. n Thick dashed lines show proton gaps which 



were used in the work of Yakovlev et al. (2004) performed 
in the framework of the "standard plus exotics" sce- 
nario. In their model proton gaps are artificially enhanced 
(that is not supported by any microscopic calculations) 
just to get a better fit of the data. We use their "lp" 
model. Neutron 3P 2 gaps presented in Fig. ^ (thick dash 
-dotted lines) are the same, as those of "3nt" model of 
Yakovlev et al. (2004) We will call this choice, the model 



I. Thin lines show ISq proton and 3P 2 neutron gaps from 
Takatsuka fc Tamagaki (2004)1 for the model AV18 by 
Wiringa et al. (1995) (we call it the model II). We take 
the same ISq neutron gap in both models I and II (thick 



solid line), as it was calculated by Ainsworth et al. (1989) 



and was previously used by Schaab et al. (1997) within 
the cooling code. [Blaschke et al. (2004) have used models 



Fig. 2. The relation between the inner crust temper- 
ature and the surface temperature for different mod- 
els. Dash-dotted curves indicate boundaries of the un- 
certainty band. Notations of lines are determined in the 



legend. For more details see Blaschke et al. (2004) and 



Yakovlev et al. (2003b) 



Recently Schwenk & Friman (2004) have argued for a 
strong suppression of the 3P 2 neutron gaps, down to val- 
ues J5 10 keV, as the consequence of the medium-induced 
spin-orbit interaction. They included important medium 
effects, as the modification of the effective interaction of 
particles at the Fermi surface owing to polarization con- 
tributions, with particular attention to spin-dependent 
forces. In addition to the standard spin-spin, tensor and 
spin-orbit forces, spin non-conserving effective interac- 
tions were induced by screening in the particle-hole chan- 
nels. Furthermore a novel long-wavelength tensor force 
was generated. The polarization contributions were com- 
puted to second order in the low-momentum interaction 
Mowk- These findings motivated Blaschke et al. (2004) 



to 



suppress values of 3P 2 gaps shown in Fig. [T] by an extra 
factor /(3P 2 ,rt) = 0.1. Further possible suppression of the 
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3P2 gap is almost not reflected on the behavior of the 
cooling curves. 




i°g(t[yr]) 



Fig. 3. Fig. 21 of Blaschke et al. (2004) Gaps are from 



Fig.nf° r model II. The original 3/2 neutron pairing gap 
is additionally suppressed by a factor /(3P2, n) = 0.1. The 
T s — 71 n relation is given by "our fit" curve of Fig.El Here 
and in all subsequent figures the value T s is the red-shifted 
temperature. NS masses are indicated in the legend. For 
more details see Blaschke et al. (2004) 



Contrary to expectations of Schwenk & Friman (2004) 
a more recent work of Khodel et al. (2004) argued that 



the 3P2 neutron pairing gap should be dramatically en- 
hanced, as the consequence of the strong softening of 
the pion propagator. According to their estimation, the 
3P2 neutron pairing gap is as large as 1 4- 10 MeV in a 
broad region of densities, see Fig. 1 of their work. Thus 



results of calculations of Schwenk fc Friman (2004)| and 
Khodel et al. (2004)[ which both had the same aim to in- 
clude medium effects in the evaluation of the 3/2 neutron 
gaps, are in a deep discrepancy with each other. 

Note that in order to apply these results to a broad 
density interval both models may need further improve- 
ments. The model of Schwenk & Friman (2004) was de- 
veloped to describe not too high densities. It does not in- 
corporate higher order nucleon-nucleon hole loops and the 
A isobar contributions and thus it may only partially in- 
clude the pion softening effect at de nsities Z no. Contrary, 
the model of |Khodel eta l. (2004) uses a simplified an- 
alytic expression for the effective pion gap ui 2 (k m ) = 
minfc[— Z? _1 (cij = 0, fc = fc m )], where D~ x is the in-medium 
pion propagator, valid near the pion condensation criti- 
cal point, if the latter occurred by a second order phase 
transition. The latter assumption means that uj 2 (k rn ) is 
assumed to be zero in the critical point of the phase tran- 
sition. Outside the vicinity of the critical point the pa- 
rameterization of the effective pion gap that was used can 
be considered only as a rough interpolation. Actually the 
phase transition is of first order and evaluations of quan- 



tum fluctuations done by Dyugaev (1982) show that the 
value of the jump of the effective pion gap in the critical 
point is not as small. Moreover repulsive correlation con- 
tributions to the NN amplitude have been disregarded. 
In the pairing channel under consideration, already out- 
side a narrow vicinity of the pion condensation critical 
point, the repulsion originating from the NN correlation 
effects may exceed the attraction originating from the pion 
softening. Notice that, if the pairing gap enhancement oc- 
curred only in a rather narrow vicinity of the pion con- 
densation critical point, it would not affect the results of 



Blaschke et al. (2004) In the latter work two possibilities 
were considered: i) a saturation of the pion softening with 
increase of the baryon density resulting in the absence 
of the pion condensation and ii) a stronger pion soften- 
ing stimulating the occurrence of the pion condensation 
for n > n c ~ 3 uq. In both cases the effective pion gap 
was assumed never approaching zero and undergoing a 
not too small jump at the critical point from a finite pos- 
itive value (ui 2 ~ 0.3 m 2 , m w is the pion mass) to a finite 
negative value (ui 2 ~ —0.1 m 2 ). The reason for such a 
strong jump is a strong coupling. If it were so, a strong 



softening assumed by Khodel et al. (2004) would not be 
realized. However, due to uncertainties in the knowledge 
of forces acting in strong interacting nuclear matter and 
a poor description of the vicinity of the phase transition 
point we can't exclude that the alternative possibility of a 
tiny jump of the pion gap exists. Therefore we will check 
how these alternative hypotheses may work within our 
"nuclear medium cooling" scenario. Thus avoiding fur- 
ther discussion of the theoretical background of the mod- 
els, in this note we investigate the possibility of a signifi- 
cantly enhanced 3P2 neutron pairing gap and of a partially 
suppressed proton ISq gap, as it has been suggested by 



Khodel et al. (2004) To proceed in the framework of our 
nuclear medium cooling" scenario we introduce the en- 
hancement factor of the original 3P2 neutron pairing gap 
/(3P2, n), and a suppression factor of the proton ISq gap 
f(lSo,p). We do not change the neutron ISq gap since in 
any case its effect on the cooling is minor. 

Generally speaking, the suppression factors of super- 
fluid processes are given by complicated integrals. As it 
was demonstrated by Scdrakian (2 005)| on the example of 
the DU process, these integrals are, actually, not reduced 
to the so called i?-factors, see Yakovlev et al. (2004) 



However, for temperatures essentially below the critical 
temperature the problem is simplified. With an exponen- 
tial accuracy the suppression of the specific heat is gov- 
erned by the factor £ n „ for neutrons, and £ pp for protons: 



~ exp[-A u (T)/T], for T < T c , 



n,p, 



(1) 



and £a = 1 for T > T c ,, T C i is the corresponding crit- 
ical temperature. We do not need a higher accuracy to 
demonstrate our result. Therefore we will use these sim- 
plified factors. 

For the emissivity of the DU process the sup- 



pression factor is given by 



SCO 



Lattimer et al. (1991) Suppression factors for 
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two nucleoli processes follow from this fact and 
from the diagrammatic representation of differ- 
ent processes within the closed diagram technique 
Voskresensky & Senatorov (1987) 



by 

Znn 



These are: 
branch of the 



min{^ ran , for the neutron 

MU process (and for the medium modified Urea process, 
MMU); £p P ■ mm{£„ n , £ P p} for the corresponding proton 
branch of the process; £ 2 n for the neutron branch of 
the (medium modified) nucleon bremsstrahlung (MnB) 
and £pp for the corresponding proton branch of the 
bremsstrahlung (MpB). Thus, for A nn 3> A pp both 
neutron and proton branches of the MMU process are 
frozen for T <C T cn due to the factors and £, pp £, nn , 
respectively. 

The resulting cooling curves depend on the T in — T s 
relation between internal and surface temperatures in 
the envelope. Fig. [21 shows uncertainties existing in this 
relation. Calculation is presented for the canonical NS: 
M = 1.4M Q , R = 10 k m with the crust model HZ90 of 
Yakovlev et al. (2003b) Below we will show that a mini- 



mal discrepancy with the data is obtained with "our fit" 
model. Using other choices like "Tsuruta law" (T 8 Tsur = 
(10Ti n ) 2 / 3 , where T s and T; n are measured in K) only in- 
creases the discrepancy. To compare results with "our fit" 
model we use the upper boundary curve, "77 = 4 • 10~ 16 " 
and the lower boundary curve "77 = 4 • 10~ 8 " . In Fig|2 
we also draw lines 77 = 1 ■ 10~ 14 and r) = 1 • 10~ n 
as they are indicated in the corresponding Fig. 2 of 



Yakovlev et al. (2003b) In reality the selection of 77 = 
4 • 10 -8 and 77 = 4 • 10 -16 as the boundaries of the 
uncertainty-band seems to be a too strong restriction, see 



Yakovlev et al. (2003b) The limit of the most massive he- 
lium layer is achieved for 77 ~ 10~ 10 . On the other hand 
the helium layer begins to affect the thermal structure only 
for 77 > 10~ 13 . Thus one could exploit 10~ 13 < 77 < 10~ 10 , 
as a Ti n — T s band. We will use a broader band, as it is 
shown in Fig[2 By this we simulate effect of maximum 
uncertainties in the knowledge of the Ti n — T s relation. 

We present Fig. 21 of Blaschke et al. (2004) now Fig. 
|2| Cooling curves shown in this figure were calculated us- 
ing "our fit" model of the crust, demonstrated by the 
solid curve in Fig. Here and in the corresponding fig- 
ures below the surface temperature is assumed to be red- 
shifted, as it is inferred by the observer from the radi- 
ation spectrum. Gaps are given by the model II of Fig. 

However, the 3P2 gap is additionally suppressed by a 
factor f(3P2,n) = 0.1, as motivated by calculations of 



Schwenk & Friman (2004) If we took the original 3P2 gap 
of the model II, we would not succeed to describe the 
data. The cooling then would be too fast, see Fig. 22 of 



Blaschke et al. (2004) Now we will check the possibility 



of ultra- high 3P2 neutron pairing gaps, as motivated by 



Khodel et al. (2004) 



In case when neutron processes are frozen the most ef- 
ficient process is the pPBF process, p — ► pvv, for T < T cp . 
Taking into account medium effects in the weak coupling 
vertex we use the same expression for the emissivity of 



this process as has been used by Voskresensky (2001) 



Blaschke et al. (2004) 



e„[(M) P PBF] ~ 10 



■ 2 „ 



m *N 


PF.p 






m N 


_PF,n(n )_ 




MeV 



7 



T 



-1 1/2 



*pp 



I'll 



erg 
cm 3 sec 



T < T a 



(2) 



We point out that this process contributes only be- 
low the critical temperature for the proton pairing. 
Inclusion of medium effects greatly enhances the ver- 
tex of this process compared to the vacuum vertex. 
Due to that a factor ~ 10 2 arises, since the pro- 
cess may occur through nn~ x and ee _1 correlation 
states (—1 symbolizes the particle hole), with subse- 
quent production of vv from the nn~ 1 vv and ee~ 1 vv 
channels rather than from a strongly suppressed chan 
nel pp~ x vv, see Voskresensky & Senatorov (1987) 



|Senatorov fc Voskresensky (1987)| |Migdal et al. (1990) 



Voskresensky et al. (1998)f 



Leinson (2000) 



[Voskresensky (2001) Relativistic corrections incor- 
porated in the description of the pp~ 1 vv vertex also 
produce an enhancement but quite not as strong as 
that arising from medium effects in nn~ x vv and ee~ x vv 
channels. We point out that we see no arguments not to 
include these corrections and we pay attention to only a 
moderate dependence of the result on the uncertainties 
in the knowledge of the strong interaction. 

We also present here an explicit expression for 
the emissivity of the proton branch of the nucleon 
bremsstrahlung including medium effects, MpB, pp — > 
ppvv. In case of suppressed neutron 3P2 gaps this pro- 
cess contributed much less than several others. However, 
in case when neutron processes are frozen, the pp — > ppvv 
process becomes the dominating process for T cn > T > 
T cp . The emissivity of the pp — > ppvv reaction takes the 



form (see Voskresensky & Senatorov (1986) for more de- 
tails) 



e(MpB) ~ 10 23 ^p/ P p 



(™kY ( — 



5/3 



v 5/3 r 2 -p4rp8 
1 P L w L s- L 9 

erg 



cm" sec 



(3) 



T 9 = T/10 9 K, m* N is the effective nucleon mass, T w ~ 1, 
and T s ~ 1/[1 + C(n/n ) 1/3 ], C ~ 1.4 H- 1.6, take into 
account TV TV correlations in weak and strong interaction 
vertices, respectively, Y p — n p /n is the proton to nucleon 
ratio. We for simplicity assumed that the value k = k m , 
at which the effective pion gap uP(k) gets the minimum, 
is rather close to the value of the neutron Fermi mo- 
m entum Pf,w (as it follows from the microscopic analysis 
of Migdal et al. (1990) I . To simplify the consideration we 
take the same value of the effective pion gap for the given 
process as that for the MMU process (although in gen- 
eral case it is not so, and thus the result (J3J) proves to be 
essentially model dependent), cf. Blaschke et al. (2004) 

K_ ( PF,n \ 5 g 
64 \PF :P J PF,n 



I pp ~ . . 



(4) 



H. Grigorian and D.N. Voskresensky: Cooling of Neutron Stars and 3P2 Neutron Gap 



5 



We have checked that for T < T cp for the pairing gaps 
under consideration the MpB reaction contributes signif- 
icantly less than the pPBF process. It could be not the 
case only in a narrow vicinity of the pion condensation 
critical point, if pion condensation occurred with only a 
tiny jump of the effective pion gap in the critical point. 
However, even in this case there are many effects which 
could mask this abnormal enhancement. 

In case of frozen neutron degrees of freedom the specific 
heat is governed by protons and electrons: 

c p ~ 10 20 K/m») (V"o) 1/3 fp P T 9 erg cm^K" 1 , (5) 



6 x 10 19 (n e /n ) 2/3 T 9 erg cm" 3 K - 



(G) 



Here, we again suppress a contribution to the specific 
heat of a narrow vicinity of the pion condensation crit- 
ical point due to the fact that in our scenario (see Fig. 
1 of Blaschke et al. (2004) I the modulus of the effective 
pion gap ui' 2 is always larger than ~ (0.1 -j- 0.3) m\. 
With such an effective pion gap the pion contribution 
to the specific heat is not too strong and can be dis- 
regarded in order to simplify the consideration. For the 
second order phase transition (either for a first order 
phase transition but with a tiny jump of \u> 2 \ in the crit- 
ical point), pion fluctuations would contribute stronger 
to the specific heat yielding a term oc T/uj, see 



Voskresensky fe Mishustin (1982)[|Migdal et al. (1990) 

In Figs. 0] and [S] we demonstrate the sensitivity of 
the results presented in Fig. [21 to the enhancement of the 
neutron 3P2 gap and to a suppression of the lSo proton 
gap, following the suggestion of Khodel et al. (2004) We 
start with the "our crust" model and the model II for 
the gaps, using however the additional enhancement fac- 
tor f(3P2,n) = 50 for the neutron 3P2 gap. Introducing 
factors f(lSo,p) = 0.1 and f(lSo,p) = 0.5 we test the 
sensitivity of the results to the variation of the 15*0 proton 
gap. We do not change the value of the lSo neutron gap 
since its variation almost does not influence on the cool- 
ing curves for NS's with masses M > 1 M , that we will 
consider. 

Comparison of Figs. 01 - El shows that in all cases NS's 
with masses M > 1.8 M cool similar in spite of the fact 
that 3P2 neutron and lSo proton gaps are varied in wide 
limits. This is because 3P~2 neutron and lS'o proton gaps 
disappear at the high densities, being achieved in the cen- 
tral regions of these very massive NS's, see Fig. ^ Thus 
these objects cool down similar to non-superfluid objects. 
Extremely rapid cooling of stars with M > 1.84 M is 
due to the DU process, being very efficient in the normal 
matter. Thereby we will notice that the cooling curves are 
very sensitive to the density dependence of the gaps. The 
difference in the cooling of NS's with M < 1.8 M in cases 
presented by Figs. 01 and is the consequence of different 
values of proton gaps used in these two calculations. This 
difference is mainly due to the pPBF processes. The larger 
the proton gap, the higher is the emissivity. 

We checked that for stars with M < 1.6 M for T < 
T cn for the 3P2 neutron pairing, a complete freezing of 



neutron degrees of freedom occurs already for /(3P2, n) > 
20. Contributions to the emissivity and to the specific heat 
involving neutrons are fully suppressed then. For heavier 
stars (M > 1.6 M ) a weak dependence on the value of 
the factor /(3P2,n) still remains even for /(3P2,?i) > 100 
but the corresponding cooling curves lie too low to allow 
for an appropriate fit of the data. This difference between 
cooling of stars with M < 1.6 M and M > (1.6-=-1.7) M Q 
is due to the mentioned density dependence of the neutron 
3P2 gap. The latter value smoothly decreases with increase 
of the density reaching zero for n > 4.5 no (the density 
4.5 no is achieved in the center of a NS of the mass M = 
1.7 Mq). At densities slightly below 4.5 no the gap is 
rather small. Therefore for stars with M > (1.6 -j- 1.7) M 
the scaling of the gap by a factor /(3P2,n) changes the 
size of the region where the gaps may affect the cooling. 
For stars with M ^ 1.6 M gaps have finite values even 
at the center of the star. Thereby there exists a critical 
value of the factor /(3P2,n), such that for higher values 
of /(3P2,n) the cooling curves are already unaffected by 
its change. 

Figs. 0] and demonstrate that we did not succeed 
to reach appropriate overall agreement with the data 
getting too rapid cooling. If we used the effective pion 
gap that allowed for the second order pion condensa- 
tion transition, as Khodel et al. (2004) assumed, we would 



get much more rapid cooling that would disagree even 
more with the data. The cooling of the old pulsars is not 
explained in all cases. Although the heating mechanism 
used by Tsuruta (2004) may partially help in this respect, 
the discrepancy between the curves and the data points 
seems to be too high, especially in Fig. We see that 
in our regime of frozen neutron processes a better fit is 
achieved in Fig. 01 i.e., for a stronger suppressed proton 
gap (for f(lSo,p) = 0.1). Actually we note that the dis- 
crepancy is even more severe, since to justify the idea of 



Khodel et al. (2004) we should exploit a softer pion prop- 



agator. Only a strong softening of the pion mode might 
be consistent with significant increase of the neutron 3P2 
gap. On the other hand such an additional softening would 
immediately result in a still more rapid cooling. The work 
oflVoskresensky et al. (2000)|discussed the possibility of a 



novel very efficient process with the emissivity e v oc T 5 
that would occur due to non-fermi liquid behavior of the 
Fermi sea in a narrow vicinity of the pion condensation 
critical point at the assumption of a strong pion softening. 
If we included this very efficient process, the disagreement 
with the data could be strongly enhanced. An enhance- 
ment of the specific heat due to pion fluctuations within 
the same vicinity region of the pion condensation point 
can't compensate the acceleration of the cooling owing to 
the enhancement of the emissivity. Khodel et al. (2004) 
used the value n c — 2 no for the critical density of the 
pion condensation. In case of the Urbana-Argonne equa- 
tion of state that we exploit here (we use the HHJ fit of 
this equation of state that removes the causality problem, 
seefBlaschkc et al. (2004) for details) the density n = 2 no 
is achieved in the central region of a NS with the mass 
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M ~ 0.8 M Q . This means that all NS's with M > 0.8 M Q 
would cool extremely fast and would not be seen in soft 
X rays. 




Fig. 4. Cooling curves according to the nuclear medium 
cooling scenario, see Fig.|3 Gaps are from Fig.^for model 
II but the 3P2 neutron pairing gap is additionally en- 
hanced by a factor /(3P2,n) = 50 and the ISo proton 
gap is suppressed by f(lSo,p) =0.1. The T s — T; n relation 
is given by "our fit" curve of Fig. [5] 



and 0.1 -j- 0.5 respectively. We verified that the variation 
of f(3P2,n) and f(lSa,p) factors in the whole mentioned 
range done within our parameterization of the effective 
pion gap does not allow to improve the picture. In all cases 
we obtain too fast cooling. To demonstrate this in Fig. H3 
we show the cooling of a 1.4 Mq star for different values 
of the /(3P2,n) factor. The factor f(lSo,p) is taken to be 
0.1. We see that for /(3P2,n) < 15 -j- 20 the curves rise 
with the increase of /(3P2,n) factor. For /(3P2,n) > 20 
the curves do not depend on /(3P2, n). 
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Fig. 5. Cooling curves according to the nuclear medium 
cooling scenario, see Fig.|3 Gaps are from Fig.^for model 
II but the 3P2 neutron pairing gap is additionally en- 
hanced by a factor /(3P2,n) = 50 and the ISq proton 
gap is suppressed by f(lSo,p) =0.5. The T s — T m relation 
is given by "our fit" curve of Fig. [5J see also Fig. 0] 

Actually, we checked the whole interval of variation 
of /(3P2,n) and f(lSo,p) factors in the range 1 -j- 100 



Fig. 6. Cooling curves of the neutron star with the mass 
1.4 Mq according to the nuclear medium cooling scenario, 
see Fig. [3J Gaps are from Fig. ^ for model II but the 3P2 
neutron pairing gap is additionally enhanced by different 
factors /(3P2,n) (shown in Figure) and the ISq proton 
gap is suppressed by f(lSo,p) — 0.1. The T s — T- m relation 
is given by "our fit" curve of Fig. 

To check how the results are sensitive to uncertain- 
ties in our knowledge of the value that determines 
the strength of the in-medium effect on the emissivity 
of the MpB process we multiplied © by a pre-factor 
/(MpB) that we varied in a range /(MpB) = 0.2 -j- 5. 
In agreement with the above discussion, for /(MpB) < 1, 
for temperatures logT s [K] > 5.9 the cooling curves are 
shifted upwards. Opposite, for /(MpB) > 1, for tempera- 
tures logP s [K] > 5.9 the cooling curves are shifted down- 
wards. However independently of the value /(MpB) for 
logT s [K] < 5.9 curves are not changed. Thus it does not 
allow to diminish the discrepancy with the data. 

Now we will check the efficiency of another choice of 
the gaps, as motivated by the model I, thick lines in Fig. 
n Compared to the model II the model I uses an artifi- 
cially enhanced proton gap. Thereby, one can expect that 
the model I is less realistic than the model II. Also we pay 
attention to a different density dependence of the pro- 
ton gap (it cuts off for densities n Z 3rto in the model 
I) compared to that given by the model II. However, as 
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we have mentioned, uncertainties in existing calculations 
of the gaps are very high. Thus it is worthwhile to check 
different possibilities. Since the mentioned parameteriza- 
tion has been used by one of the groups working on the 
problem of cooling of NS's, see Yakovlev et al. (2004) we 



will consider consequences of this possibility as well. Fig.[7| 
demonstrates our previous fit of the data within the model 
I, but for the original 3P 2 neutron gap being suppressed 
by /(3P 2 ,n) = 0.1 (see |Blaschke et al. (2004) I. 




Fig. 7. Fig. 15 of |Blaschke et al. (2004)] Gaps from the 
model I. The original 3P 2 neutron gap is suppressed 
by /(3P 2 ,n) = 0.1. The T s - T in relation is given 
by "our fit" curve of Fig. [21 Other notations see in 
Blaschke et al. (2004) 



Figs. |H1 and show that within the variation of the 
gaps of the model I the discrepancy with the data is still 
stronger compared to that for the above calculation based 
on the use of the model II. 

The difference between curves shown in Figs. [H] and ED 
is less pronounced than for those curves demonstrated in 
Figs. 01 and [5j Indeed, as we have mentioned, the density 
dependence of the proton gap is different in models I and 
II. In the model II the proton gap reaches up to higher 
densities than in the model I (in the latter case the gap 
is cut already for n ^ 3no). Thus in the case shown by 
Figs. 03 and [5] a non-superfluid core begins to contribute 
already for smaller values of the star mass. Therefore in 
both figures the corresponding cooling curves are almost 
the same for M > 1.6 A/©. 

The dependence of the results on the different choices 
of the T s — Ti n relation is demonstrated by Figs. llOlllll for 
gaps based on a modification of the model II and by Figs. 
ElEDf 01 ' gaps based on a modification of the model I. For 
this demonstration we first took the upper boundary curve 
77 = 4- 10~ 16 in Fig. 13 We show that these choices however 
do not allow to improve the fit. Comparing Figs. 1 101 and l4l 
(and II Hand 15)) based on the very same modification of the 




Fig. 8. Cooling curves according to the nuclear medium 
cooling scenario, see Fig.0 Gaps are from Fig.^for model 
I but the 3P 2 neutron pairing gap is additionally enhanced 
by a factor /(3P 2 ,n) = 50 and the ISq proton gap is 
suppressed by f(lSo,p) = 0.1. The T s — T; n relation is 
given by "our fit" curve of Fig. [21 




Fig. 9. Same as Fig. [SI but for the original ISq proton gap 
suppressed by f(lSo,p) = 0.5. 



model II we see that with the "our fit" crust model the 
deviation from the data points is less pronounced. Basing 
on the model I we compare Figs. I12landl81 (and 1131 and !§)) 
and we arrive at the very same conclusion. 

In Fig. 1141 we use the lower boundary curve r\ = 
4.0 • 10~ 8 of the Fig. [21 We further demonstrate that 
the selection of a different choice of the T s — T- m relation 
within the band shown in Fig. [21 does not allow to dimin- 
ish discrepancy with the data. Contrary, this discrepancy 
just increases compared to that demonstrated by "our fit" 
model. Indeed, the cooling evolution for times t ^ 10 5 yr 
(T s ^ 10 6 K) is governed by neutrino processes. Thus the 





Fig. 10. Same as in Fig. 0] but for the crust model r\ — Fig. 12. Same as in FiglHl but for the crust model r\ 
4.0 ■ 1(T 16 . 4.0 ■ 1(T 16 . 





log(t[yr]) 



log(t[yr]) 



Fig. 11. Same as in FigEU but for the crust model r\ = Fi g- 13- Same as in Fig. for the crust model rj = 4.0 
4.0 -10- 16 . 10^ 16 - 



higher Tj n , the larger T s is. The slowest cooling is then ob- 
tained, if one uses the lower boundary curve r\ = 4.0 • 10 -8 
of Fig. [21 The evolution of NS's for times t 10 5 yr begins 
to be controlled by the photon processes. In the photon 
epoch (t 3> 10 5 yr) the smaller the T s value, the less ef- 
ficient the radiation is. Thus for t 3> 10 5 yr the slowest 
cooling is obtained, if one uses the upper boundary curve 
77 = 4.0 • 10~ 16 of Fig. |2 The "our crust" curve just sim- 
ulates the transition from the one limiting curve to the 
other demonstrating the slowest cooling in the whole tem- 
perature interval shown in the figures. 

We point out that in all cases the data are not ex- 
plained within the assumption of an enhanced neutron 
3P2 gap (for f(3P2,n) > 1) and of a partially suppressed 
ISo proton gap (for f(lS ,p) = 0.1 0.5). 



Our aim was to consider the possibility of large ZP2 
gaps within the same "Nuclear medium cooling" sce- 
nario of Blaschke et al. (2004) that well described the 



cooling data in opposite assumption of suppressed 3P2 
gaps. Therefore in the present work we did not incor- 
porate possibilities of internal heating for old pulsars 
(see Tsuruta (2004) I and of existence of quark cores 



(see Blaschke et al. (2001) Grigorian et al. (2004) and 
refs therein). 

The main problem with the given scenario is the fol- 
lowing. At the frozen neutron contribution to the specific 
heat and to the emissivity the region of surface tempera- 
tures T s > 10 6 K is determined by proton processes. The 
most efficient among them is the pPBF process. For proton 
gaps, which we deal with, the pPBF process proves to be 
too efficient yielding too rapid cooling. Thus at least sev- 
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log(t[yr]) 

Fig. 14. Same as in Fig. 0] but for the crust model rj = 
4.0- 1(T 8 . 

eral slow cooling data points (at least data for old pulsars) 
are not explained. Note that some works ignore the men- 
tioned above medium induced enhancement of the pPBF 
emissivity that results in 10-100 times suppression of the 
rate. We omitted this possibility as not physical one. 
Concluding, we have shown that the "Nuclear medium 



cooling" scenario of Blaschke et al. (2004) fails to appro- 
priately fit the NS cooling data at the assumption of a 
strong enhancement of the 3P2 neutron gaps (we checked 
the range /(3P2,n) = 1 100) and for moderately sup- 
pressed ISq proton gaps (for f(lSo,p) = 0.1-^0.5). On the 
other hand the very same scenario allowed us to appropri- 
ately fit the whole set of data at the assumption of a signif- 
icantly suppressed 3P~2 neutron gap (for f(3P2,n) ~ 0.1). 
We observed an essential dependence of the results not 
only on the values of the gaps but also on their density 
dependence. We used the density dependence of the gaps 
according to the models I and II. The latter model is sup- 
ported by microscopic calculations. We excluded an at- 
tempt to artificially fit the density dependence of the gaps 
trying to match cooling curves with the data. Although 
such an attempt could improve the fit, this way seems us 
rather not physical and we did not follow it. However we 
will greatly encourage further attempts of microscopic cal- 
culations of the gaps, which would take into account most 
important medium effects. With carefully treated gaps one 
could return to the simulation of the NS cooling. 
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